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Abstract. - We reanalyze the data on the in-plane far infrared conductivity of underdoped or- 
tholl YBa2Cu3O6.50 (YBCOg.so) in terms of a model in which a pseudogap opens on part of the 
Fermi surface with the remaining ungaped piece proportional to the temperature. The motivation 
for our model comes from recent angle-resolved photoemission spectroscopy data in underdoped 
Bi2Sr2CaCu208+i (Bi-2212), which have revealed the existence of arcs on the Fermi surface. We 
find the optical data to be consistent with arc formation. In addition we find some evidence 
that the electronic states lost below the pseudogap energy A pg are recovered in the energy region 
immediately above it at least for temperatures near T c . 
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In going from the overdoped side of the cuprate phase 
diagram to the underdoped side, a new phenomenon 
known generally as the formation of a pseudogap is seen 
[1]. The pseudogap manifests itself in many properties and 
generally has the signature of a decrease in the electronic 
density of state (DOS) around the Fermi energy. There is 
no consensus as to its exact microscopic origin. Some ob- 
servations are consistent with a preformed pair model [2] 
with the pseudogap associated with the binding energy of 
\ the pairs at the pseudogap temperature T* and the super- 
conducting temperature T c the point at which phases lock 
in and long range order results. Another influential model 
is a possible competing order such as d-density [3-8] wave 
which set in at T* and coexists with the superconducting 
order parameter below T c . In this model it may be rea- 
sonable to expect two distinct gap scales and some recent 
data [9, 10] have pointed to this possibility. 

In the absence of consensus on microscopic origin, phe- 
nomenological models have played an important role in 
understanding and correlating data. An example is the 
specific heat data [11, 12]. The observed trend from over- 
doped to underdoped can be understood [11,12] over sev- 
eral families of cuprates by simply introducing a gap A pg 
in the DOS which obeys a mean field temperature depen- 
dence. More recently, analysis of angle-resolve photoe- 
mission spectroscopy (ARPES) data have given a more 
detailed picture [13] of pseudogap formation on the Fermi 



surface as a function of temperature. A finite pseudogap 
first opens at T* but forms only in a region around the 
antinodal direction leaving an ungaped Fermi surface arc 
around the nodal direction of length proportional to the 
reduced temperature t = T/T*. Here we confine our at- 
tention to the normal state region (above T c ) and study 
how the formation of Fermi surface arcs modifies optical 
properties. 

Optical conductivity measurements have given a wealth 
of information on electron dynamics in the cuprates. 
While early on the pseudogap was seen clearly as a real 
gap in the real part of the c-axis [14, 15] optical conduc- 
tivity, its signature in the ab-plane data was less clear and 
largely confined to a general observation that the opti- 
cal scattering rate seem to show an additional decrease 
with decreasing temperature at low frequency as com- 
pared with overdoped samples [16]. Only recently has 
the situation been further clarified when the real part of 
the optical self-energy was analyzed in more detail [17]. 
Underdoped samples show a characteristic evolution with 
decreasing temperature quite distinct from the overdoped 
case [18-20]. A hat like peak structure develops at an 
energy of the order of the pseudogap energy, clearly seen 
above a smooth structureless background. It was recog- 
nized that this easily identified [17] feature could be un- 
derstood in a simple DOS model with a gap below A pg 
and the missing states piled up in the energy region just 
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above the pseudogap which we will refer to here as the 
recovery region. These modification of the DOS show up 
very directly in the real and imaginary part of the optical 
self-energy [21-24]. 

Consider coupling to a single Einstein mode of energy 
we- Without a pseudogap (constant electronic DOS) the 
optical scattering rate 1/t° p (ui) will be zero till w = uie 
after which it rises from zero according to the law (w — 
we)/w [25]. With a pseudogap it will be zero till ui = 
u>e + A pff and rise after this according to (w — uje — jw 
[20]. On the other hand if a recovery region is included 
above uj — A pg to conserve electronic states, there will 
be additional scattering in the recovery region due to the 
extra available states. This will make 1/t° p (uj) rise more 
rapidly than it would have without the extra states. This 
rapid rise in \/t° p (uj) translates into the formation of a 
logarithmic like structure in Kramers-Kronig transform. 
Monitoring the growth of this structure in the real part 
of the optical self-energy Yf^(uj) gives us a means to trace 
pseudogap formation as does the steepness of 1/t° p (uj) 
above the energy of the gap plus boson. 

There have been many studies of the boson spectra of 
the cuprates as revealed by tunneling [26], ARPES [27] 
and optics [20, 28-33] in which a spectral density func- 
tion I 2 x(w) is recovered, where / is a coupling constant 
between electrons and a boson and is the bosonic 

spectral density. For optimally and overdoped cases, cou- 
pling to the sharp resonance at an energy corresponding 
to the spin-1 collective mode measured in neutron scatter- 
ing is seen. This resonance starts to exist above T c in the 
undcrdoped case and its intensity (area under the peak) 
increases with decreasing temperature. For the specific 
case of ortholl YBCO6.50 this has been studied by Hwang 
et al. [20]. 

In this paper we reanalyze the data on ortholl 
YBCO6.50 allowing for the opening of a pseudogap which 
covers more of the Fermi surface, as the temperature is 
lowered below T* . We ask whether or not the data is con- 
sistent with a linear decrease in the length of the remaining 
arcs as a function of the reduced temperature t — T/T* as 
it is seen in the ARPES data in underdoped Bi-2212 [13]. 
Our analysis of the optical data proceeds through the op- 
tical scattering rate 1/t° p (T, uj). For a system with an en- 
ergy dependent self consistent electronic density of state 
N(u>) which is related to the imaginary part of the Green's 
function G(k,uj) by N(w) = l/^^2 k [— ImG(k, w)] the op- 
tical scattering rate can be written approximately [22] as 

n POC pOO 

1/t° p (X» =— / dm 2 X (^) deN(e-fl)x 

u Jo J-00 (1) 

[n(n)+f(n-e)][f(e-w)-f(e + w)] 

Here n(Q) = l/[e^ n — 1] is the Bose-Einstein thermal 
distribution, /(e) = l/[e p(ji ~ e ^ + 1] the Fermi-Dirac with 
j3 the inverse temperature and fi the chemical potential, 
and I 2 x{w) the electron-boson spectral density. In Eq. 
1 N(e) is the symmetrized electronic density of states 
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Fig. 1: (color online) Top frame, optical scattering rate 
1/t° p (lu) in cm -1 as a function of u> also in cm -1 for the 
eight temperatures shown (normal state). The heavy solid 
squares are fits to the data (solid lines). The inset shows the 
frequency dependent of N(w,t) for t — 0.5, with quadratic 
frequency dependence below A pg and a constant above it for 
ui 6 (A P9 , 2A ra ). Electronic states are conserved. Bottom 
frame, same as top but previous fit to the data obtained in 
reference [20] without including a recovery region above to 
pseudogap energy which is chosen to conserve state. Middle 
frame, difference between top frame and scattering rates ob- 
tained when recovery region is excluded, keeping other param- 
eter the same (as top frame fits). 

[N(e) + N(-e)]/2. Eq. 1 reduces properly to that of 
Shulga et al. [34] in the limit of a constant density of 
states. It further corresponds to the formula of Allen [35] 
at zero temperature. The data for 1/t° p (T, uj) at eight 
temperatures in the normal state of ortholl YBCO6.50 are 
reproduced [20] in Fig. 1 (top and bottom frames) as 
(dark) solid lines. In all cases the data is fit (solid dia- 
monds of different colors as indicated in the figure) with 
an electron-boson spectral density I 2 x(w) consisting of 
a broad Millis-Monien-Pines (MMP) background [36] of 
the form I uj/(uj 2 + uj 2 j), where I is a coupling between 
spin fluctuation and charge carriers and uj s j is a typical 
spin fluctuation energy. In addition a sharp resonance 
peak modelled by I p /[V2^(d/2.35)}e- ( - bJ - UJ ^ 7[2(d/2. 35 ) 2 ] 
(where I p is the area under the peak, ui v is its center 
frequency, and d is its width) is included and its po- 
sition in energy fixed at 31 meV chosen to agree with 
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Fig. 2: (color online) Model electron-boson spectral density 
I 2 x(u) at various temperatures, obtained from our fits to the 
scattering rate data assuming an optical resonace peak at 31 
meV superimposed on an MMP background. 



the energy of the measured spin-1 resonance of inelas- 
tic neutron scattering [37]. In addition a pseudogap is 
included with A pg = 350 cm -1 fixed from our previous 
analysis [20] and assumed temperature independent. The 
length of the remaining ungaped arc is modelled through 
the self consistent DOS N(u>) which we take to be equal to 
N(0) + {1-N(0))(uj/A ps ) 2 forw < A pg and 2(1 - A>(0))/3 
for u) £ (Apg, 2 Apg) and 1 beyond this. The value of N(0) 
is temperature dependent and assumed to be proportional 
to the reduced temperature t = T/T*. This agrees with 
the idea that the arc closes linearly in reduced tempera- 
ture as T is lowered below the pseudogap temperature T* . 



Further the lost states below uj = A pg due to the opening 
of the pseudogap are taken to reappear above A pg in the 
interval up to 2A pg . The fitting parameters to the data 
are area under the resonance peak (I p ) and the background 

The fits obtained are seen in the top frame of Fig. 1 
and are seen to be excellent. They are better than the 
earlier fits obtained in ref. [20] where a pseudogap was 
included but without arcs and recovery region. These fits 
are reproduced in the lower frame and are clearly not as 
good. In particular, in the region around ui — 1000 cm -1 
the new fits capture the narrow range of variation with 
temperature observed in the data. This effect is a di- 
rect consequence of the temperature dependence assumed 
for the DOS which itself reflects the temperature depen- 
dence of the Fermi surface arcs. When the arcs become 
smaller with reducing normalized temperature t — T/T* 
more states are lost in the pseudogap region as a results 
of a reduction in the value of N(u>). But in our model for 



N(u>) which is shown in the inset of the upper frame of 
Fig. 1, it is assumed that the lost states are to be found 
in the region between A pg and 2A pg so as to respect con- 
servation of states in the electronic DOS. The existence of 
this recovery region directly leads to more scattering i.e. 
an increase in l/r op (w) in this energy range than would 
be the case if no recovery occurred i.e. N(u>) was equal 
to its background value of 1. This effect pushes up the 
1/t° p (lu) curves more at low temperatures than at high 
temperature where the amount of extra states due to re- 
covery is reduced and finally vanishes in our model at the 
pseudogap temperature T — T* . It is therefore through 
the effect on l/r op (T, uj) of the recovery region that we 
see in optics the arcs and their temperature dependence. 
The arcs have their greatest effects at low temperature 
as can be seen in the middle frame of Fig. 1 where we 
show the contribution to the optical scattering rate due to 
the recovery states in N(uj) above A pg . At high temper- 
atures this contribution vanishes as N(u>) becomes equal 
to 1 everywhere. We note that, with increasing T away 
from T c , the effect of the recovery region in l/r op (T, uj) 
becomes small and we cease to be able to extract from the 
data any quantitative information on the arcs and on the 
recovery region above it. 

In Fig. 2 we show the results of our fits for the electron- 
boson spectral function I 2 x(uj) as a function of frequency 
and temperature. At each temperature the background 
was allowed to vary as was the area under the resonance 
peak fixed at an energy of 31 meV [37]. We see that the 
peak disappears with increasing temperature and at 250 
K and above only the background remains. Relative to 
the peak, the temperature dependence of the background 
is small. 

In Fig. 3 we show results for the area under the peak 
as a function of temperature above the superconducting 
T c = 59 K i.e. in the normal state. The points are denoted 
by (blue) solid diamond, they are also compared with our 
previous results from optics (green)open hexagon [20] and 
with neutron results (purple) solid pentagon [37] . It is seen 
that the fit to the neutron data has improved as compared 
to our previous optically derived estimates. 

Data on the optical scattering rate in ortholl YBCO6.50 
show clear evidence for the formation of a pseudogap 
which opens at T* (the pseudogap temperature) . We have 
used a model of the pseudogap with Fermi arcs with linear 
variation in reduced temperature t = T/T* as suggested 
in ARPES experiments on underdoped Bi-2212. The frac- 
tion of the Fermi surface which remains ungaped provides 
a finite but reduced value for the self consistent density of 
state at u> —0 with N(T,uj = 0) linear in t. A quadratic 
dependence of N(T, uj) in uj is assumed till uj = A pg af- 
ter which point we have assumed there is a region of in- 
creased DOS which contains all the missing states so that 
by uj — 2A pg states are conserved. This is a reasonable 
assumption and certainly supported by our fit for T near 
T c . At higher temperatures however, while the data re- 
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Fig. 3: (color online) Comparison of the area under the mag- 
netic resonance mode seen in inelastic neutron scattering ex- 
periments (solid purple pentagon) [37] and the area under the 
resonance peak at 31 meV of Fig. 2 (sold blue diamond). Also 
shown are the previous fits of ref. [20] (empty green hexagon) 
which did not include a recovery region in the electronic density 
of state model. 



mains consistent with the assumed model, details such as 
the extent in energy and exact amount of recovery could 
not be quantitatively pinned down. With this model we 
calculate the optical scattering time t° p (lu) at several tem- 
peratures in the normal state with electron-boson spectral 
density consisting of an MMP background and resonance 
peak set at 31 meV to agree with the available inelastic 
neutron scattering data on the local magnetic susceptibil- 
ity [37]. Two parameters, the amplitude of the background 
and the area under the resonance peak are varied to get 
a best least square fit to the data. The quality of the fit 
obtained in this way is considerably superior to an earlier 
fit also based on a pseudogap model but without recov- 
ery and Fermi surface arcs. We conclude that the data is 
more consistent with arcs than without. The effect of the 
arc shows up clearly in the region just above the main rise 
in scattering rate where all the curve bunch up together. 
This is properly captured in our fits while it is missing in 
the previous analysis. Also our new fits agree better with 
the neutron studies as to the temperature dependence of 
the area under the resonance peak which is one of the 
important parameter varied in our least square fit. 
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